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Recently, Weniger (delta sequence) method has been proposed by the au- 
thors of Ref. for resummation of truncated perturbation series in quantum 
0^ . field theories. Those authors presented numerical evidence suggesting that 

this method works better than Pade approximants when we resum a function 
with singularities in the Borel plane but not on the positive axis. We present 
here numerical evidence suggesting that in such cases the combined method 
' of Borel-Pade works better than its analog Borel- Weniger, and that it may 
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work better or comparably well in some of the cases when there are singular- 



I. INTRODUCTION 



In this letter we want to present some numerical results which allow us to compare the 
efficiency of the Borel-Pade method with that of the Borel- Weniger method for resumma- 
tion of truncated perturbation series (TPS) in some physically significant scenarios. The 
scenarios we are referring to are those when the function, which we want to find through a 
resummation, is known to have certain singularity structure in the Borel plane. If there are 
singularities on the positive axis of the Borel plane, then we implicitly assume that in such 
cases we either know the correct prescription for integration in the Laplace-Borel integral, 
or we simply adhere to a certain adopted prescription. 

(1) We will first illustrate the efficiency of the two methods on the QED example of the 
Euler-Heisenberg Lagrangian density, i.e., the one-loop fermion-induced effective action 
density in a strong uniform electromagnetic field @- [§. In this case, the solution is known, 
and its real part can be written in the following form: 
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Here, <5£ is the actual Lagrangian density induced by the one-loop fluctuations of the 
fermions in the field; g is the field-to-fermion coupling parameter (in QED it is the positron 
charge eo); m is the mass of the fermion (electron); a and b are Lorentz-invariant expressions 
characterizing the electric and the magnetic fields S and B, respectively 
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Expression ([!]) can be obtained, for example, directly by integrating out the fermionic de- 
greees of freedom in the path integral expression of the full effective action, then employing 
the proper-time integral representation for the difference of logarithms, evaluating the traces 
in the integrand, and subsequently performing Wick rotation by — 7r/4 in the plane of the 
proper-time s: ags t— > — iw + e'. We refer to for more details on the latter point. The 
perturbative expansion of the full solution (ID), in powers of a, is 
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In the case of the pure magnetic field (p.m.f.), the corresponding expressions are simpler 
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We can now use (|])-(|5|), and (0), as a laboratory for resummation methods, since the full 
(resummed) solutions ([!]) and (||) are known. Since ([!]) and @ are Laplace-Borel integrals, it 
is natural to use these examples for testing combined resummation techniques which involve 
Borel transformation. Borel transform Bl of series (f§) is 



B L (w;p) = Ci(p)w + c 3 (p)w 3 + c 5 (p)w 5 + 



(8) 



and analogously for (0). In Ref. 0, we used Borel-Pade technique for resummation, i.e., we 
applied various Pade approximants [N/M]b(w;p) to (||)[] and then employed the Laplace- 
Borel integral to obtain the resummed value 



1 [N/M]b(w;p), being ratio of polynomials in w of powers iV and M, respectively [1C], is based 
solely on the truncated perturbation series (TPS) of @ involving only terms with c n : n < N+M. 
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The integration over poles in was carried out according to the Cauchy principal value 
prescription, since the full solution ([]]) requires it.0 

Recently, the authors of |L| proposed the use of Weniger (delta sequence) transformations 
as an alternative to the use of Pade approximants, for direct resummation of truncated 
perturbation series. For a truncated perturbation series (TPS) of the form 
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where (C+j)n-i = r(C+J+ 72— l)/r(C+i) are the Pochhammer symbols and = 1 is usually 
taken. The approximant (|TTD is a ratio of two polynomials in z of power n each, and when 
expanded back in powers of z it reproduces all the terms of i^n+i]- 

The authors [I]] applied ( |TTD directly to the TPS's of 5C pcrt (a;p)/a 2 of (f|), and when re- 
expanding the approximant in powers of a they were able to predict the next coefficient in the 
series with a better precision than the one provided by the corresponding diagonal (or almost 
diagonal) Pade approximant. Further, in the case of the pure magnetic field they showed 
that the method (0), when applied directly to the TPS's in b of the induced Lagrangian 
density,[] gave better results of resummation than the corresponding Pade approximants. 

We now combine the method ( |Tl~D with the Borel transformation (f|) i— > (|8|), and compare 
the results of resummation obtained in this way with the results of the corresponding Borel- 
Pade approximants of Ref . . Formula ( |TTD is applied to the Borel transform (|]) divided by 
w. We identify z = w 2 (we thank the authors of ]15| for pointing out that this clarification was 
missing in the original version of the preprint). In the ensuing Borel- Weniger approximant, 
we integrate in the Laplace-Borel integral over the poles of the integrand with the Cauchy 
principal value prescription, just as in Borel-Pade approximant ([]), in accordance with the 
full known solution (JJ). 

The results of these calculations are presented in Figs. |T|(a)-(d), as functions of the 
electric field strength parameter a, for various values of p = b/a = 0.,0.5, 1.5,5.0. In Fig. |2| 
we present the analogous results for the case of the pure magnetic field (p.m.f.), as function of 



2 Various QCD and QED applications of the Borel-Pade approach with the principal value pre- 
scription can be found in [0]- |l3[ . The novel method of Ref. is, in addition, well suited for 
obtaining the imaginary part of 5C. 

3 The approximants (|TT|) applied to the TPS's of the series @ divided by b 2 . 
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the magnetic field parameter b. N3 and [3/4] denote the Borel-Weniger and the Borel-Pade 
resummations based on the truncated Borel transform (^) with the first four nonzero terms 
(i.e., three terms beyond the leading order); N5 and [5/6] are based on the first six terms in 
(§). Comparison with the exact solutions, also present in the Figures, shows that Borel-Pade 
is better than the corresponding Borel-Weniger, except in the case of p = 5.0 (electric field 
combined with a much stronger magnetic field). Fig. suggests that Borel-Pade is better 
than Borel-Weniger for resummation of functions whose Borel transforms have singularities 
only outside the positive axis. Further, comparison of Fig. 2 with the results of Table I of 
Ref. [Q] suggests that Borel-Pade and Borel-Weniger methods are much more efficient than 
Weniger method in resumming series with singularities in the Borel plane. Weniger method 
in the p.m.f. case is better than Pade method [|T|. 

We can also do analogous calculations for the induced energy densities SU 
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where SU = 8TT 2 5U/(m 4: a 2 ). In that case, the simple Borel transform has a double-pole 
structure on the positive real axis, and the Pade and Weniger approximants have trouble 
simulating such multiple poles adequately. Therefore, we employ a slightly modified Borel 
transform in the case of the induced energy densities 



w 2 w 4 
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which has no multiple-pole structure - all the poles are simple. The (modified) Laplace- 
Borel integral in this case is 
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where again the Cauchy principal value has to be taken, once MB u (w;p) is replaced in (|Tj 
by its Pade or Weniger approximants. For details, we refer to Ref. where Borel-Pade was 
employed also for the induced energy densities. Weniger formula ( p~l| ) is now applied to the 
modified Borel transform ( |I~6"D divided by w 2 . The results are presented in Figs. [3](a)-(d), as 
functions of a at fixed p = 0.,0.5, 1.5,5.0, respectively^ We present the solutions of Borel- 
Weniger and Borel-Pade based on the first four (N3, [4/4]) and six (N5, [6/6]) nonzero 
terms of the modified Borel transform of the energy density. We see that for the induced 



4 In the case of the pure magnetic field, the energy density is the same as the Lagrangian density, 
except for the sign change. 
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energy density the situation is less clear. In the cases p = 0, 0.5 and 5.0 the Borel-Pade 
and Borel-Weniger resummations are apparently of comparable quality, while at p — 1.5 the 
Borel-Pade appears to work better. 

We can see these trends also if we compare the perturbation coefficients predicted by 
these two methods with the exact ones. These results are written in Table | for the case of the 
Lagrangian density (predicted eg and C13) and in Table |lj for the case of the energy density 
(predicted dg and c/13) . Predictions of Borel-Pade and Borel-Weniger are of comparable 
quality in the cases of p = 0., 0.5, 5.0 for energy density and in the case of p = 5.0 for 
Lagrangian density. In other cases, predictions of Borel-Pade are better. In fact, in the 



approximant 


P = 


0.0 and p.m.f 




p = 0.5 






p = 1.5 




V 


= 5.0 


N3 


eg -- 


= 2.1666- 10" 6 


eg = 


= 3.524-10" 


-6 


Cg 


- 4.320-10" 


-4 


eg = 


596.91 


[3/4] 


eg : 


= 2.1637- 10" 6 


eg = 


= 3.648-10" 


-6 


eg 


= 5.866-10" 


-4 


eg = 


595.28 


exact 


eg = 


= 2.1644- 10" 6 


eg = 


= 3.711-10" 


-6 


eg 


= 6.166-10" 


-4 


eg = 


596.24 


N5 


C13 


= 2.2212-10" 8 


C13 


= 3.725-10" 


-8 


C13 


= 2.460-10- 


-5 


C13 = 


3823.65 


[5/6] 


C13 


= 2.2215-10" 8 


Cl3 


= 3.804-10' 


-8 


Cl3 


= 3.157-10 


-5 


Cl3 = 


3824.42 


exact 


C13 


= 2.2215-10" 8 


C13 


= 3.805-10' 


-8 


C13 


= 3.161-10- 


-5 


C13 = 


3824.45 



TABLE I. Coefficients eg and C13 of the perturbation series for the induced Lagrangian density, as 
predicted by various Borel-Weniger and Borel-Pade approximants. We include exact values for comparison. 
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TABLE II. Coefficients dg and di3 of the perturbation series for the induced energy density, as predicted 
by various Borel-Weniger and Borel-Pade approximants. For comparison, exact values are included as well. 

case p = 5.0 of the energy density, the modified Borel-Weniger is slightly, but discernibly, 
better than the modified Borel-Pade. Comparing predictions of Table | (for p = 0.0) with 
those of Tables II and III of Ref. |L] suggests strongly that the discussed Borel-Pade and 
Borel-Weniger methods are better than Weniger method in predicting the coefficients c n . 
Weniger method is better than Pade method in predicting c n 's 

(2) The second example to compare the efficiency of the Borel-Pade and Borel-Weniger 
methods will be taken from QCD, and it will have to do with the "fixing" of a pole of a 
Borel transform rather than with a resummation. We look at the Bjorken polarized sum rule 
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(BjPSR), which involves the isotriplet combination of the first moments over xej of proton 
and neutron polarized structure functions 

jf 1 dx B3 [g{ P \x B] ;Q 2 ph ) - <&)] = 1\ 9 a\ [l - S(Ql h )} . (18) 

Here, p 2 = — Q p ^< is 7* momentum transfer. At Q 2 h = 3GeV 2 where three quarks are 
assumed active (nf — 3), and if taking MS scheme and renormalization scale (RScl) Ql = Q 2 h , 
we have the following TPS of the BjPSR observable S(Q 2 h ) available |I6|| - fT7fl : 

S[2](Q 2 ph ; Ql = Q 2 ph ; cf s , cf s ) = a„(l + 3.583a + 20.215a 2 ) , (19) 
with : a = a(ln Q 2 ; cf 3 , cf s , . . .) , n f = 3 , cf s = 4.471, cf s = 20.99 . (20) 

Here we denoted by a the strong coupling parameter a = a s /ir. 

It is known from |JXS|] — ||T9| that the Borel transform B$(z) of 5* has the lowest positive pole 
at Zpoic = 1/A) = 4/9 (leading infrared renormalon) and that this pole has a much stronger 
residuum than the highest negative pole at z po \ c = — l//3 (leading ultraviolet renormalon). 
The question we raise here is: How well can Pade and Weniger approximants to the Borel 
transform Bs{z) determine the next coefficient r 3 of the term r 3 ag in the TPS (^), via the 
requirement that z po \ e = 4/9? For that, we have to know well the actual r 3 . That term 



can be determined reasonably well on the basis of two approximants discussed in pCj - the 
effective charge approximant (ECH) *4.£ ECH ^(c 3 ) with c 3 ~20., and another, also RScl- and 

1 /2 

scheme-independent approximant AJ 2 (c 3 ) with c 3 ~15.5. These two approximants give the 
correct location of the leading infrared renormalon pole, and when we expand them back in 
powers of we obtain r 3 « 129.4 and r 3 « 130.8, respectively. Therefore, we can estimate 
with high confidence the actual r 3 : r 3 = 130. ± 1. 

It is important to consider the RScl- and scheme-invariant Borel transform when we 
want to apply Pade or Weniger approximants to it, so that the predicted values of r 3 will 
be independent of the RScl- and scheme in which we work at the intermediate stage. Such 
a Borel transform has been used in [ETf] , and we use its variant Bs(z) as specified in [cf. 
Eqs. (18)-(20) there]. Such a Borel transform reduces (up to a z-dependent nonsingular 
factor) to the usual Borel transform in the approximation of the one-loop evolution. The 
resulting power expansion of Bs{z) up to ~2 3 will depend on the coefficient r 3 

B s {z) = 1 + §(7-l)y + 0.02078.. .y 2 + JL(_ 2 1.88... +jU)y 3 + 0(y 4 ) , (21) 

where 7 = 0.577... is Euler constant, and y = 2[3qz. If we apply [2/1] and [1/2] Pade 
approximants to the TPS ( |2"T| ) and demand 2; po i e = l//3 (y po \ e = 2), we obtain predictions 
r 3 = 137.0 and r 3 = 128.0, respectively. The prediction of [1/2] is significantly better, and 
this could possibly be explained with the more involved denominator structure of [1/2] in 
comparison to [2/1]. When applying to (|2"TD Weniger formula ( |TT]) with ( = 1), we 
obtain r 3 = 135.3. This is further away from the actual value of 130. ± 1. than the prediction 
of [1/2]. In both [1/2] and 82 , the denominators are polynomials of quadratic degree in z. 

To summarize this QCD example: We applied Pade and Weniger approximants to a 
(TPS of a) Borel transform of the Bjorken polarized sum rule and demanded that the 
leading infrared renormalon pole be reproduced correctly. Weniger approximant 8^ then 
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apparently gives a somewhat worse prediction for the next coefficient than the corresponding 
Pade approximant [1/2]. 
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FIG. 1. Borcl-Pade approximants ([3/4], [5/6]) and the corresponding Borel-Weniger approximants 
(N3, N5) to the induced dispersive Lagrangian density (|l|), as functions of a, for various values of p = b/a: 
(a) p = 0.0; (b) p = 0.5; (c) p = 1.5; (d) p = 5.0. The numerically exact curves are included for comparison. 



9 




C 0.2 
CD 

73 0.18 
&0.16 



FIG. 2. Borel-Pade approximants ([3/4], [5/6]) and Borcl-Wcniger approximants (N3, N5) to the 
induced dispersive Lagrangian density (^|), as functions of b, for the pure magnetic field case (5 = 0). The 
numerically exact curve is included for comparison. 
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FIG. 3. Modified Borel-Pade ([4/4], [6/6]) and the corresponding modified Borel-Weniger (N3, N5) 
approximants to the induced energy densities (|l3|), as functions of a, at fixed values of p = b/a: (a) p = 0.0; 
(b) p = 0.5; (c) p=l,5; (d) p=5.0. 
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